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Colossal Effects in Transition Metal Oxides 
can be caused by Intrinsic Inhomogeneities 
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The influence of quencfied disorder on tfie competition between ordered states separated by a 
first-order transition is investigated. A pliase diagram with features resembling quantum-critical 
behavior is observed, even using classical models. The low-temperature paramagnetic regime consists 
of coexisting ordered clusters, with randomnly oriented order parameters. Extended to manganites, 
this state is argued to have a colossal magnetoresistance effect. A scale T* for cluster formation is 
discussed. This is the analog of the Griffiths temperature, but for the case of two competing orders, 
producing a strong susceptibility to external fields. Cuprates may have similar features, compatible 
with the large proximity effect of the very underdoped regime. 

PACS numbers: 71.10.-w, 75.10.-b, 75.30.Kz 



Complex phenomena such as "colossal" magnetoresis- 
tance (CMR) in manganites and high temperature su- 
perconductivity (HTS) in cuprates have challenged our 
understanding of correlated electrons [|l| . Recent devel- 
opments unveiled a previously mostly ignored aspect of 
doped transition- metal-oxides (TMO): these systems are 
intrinsically inhomogeneous, even in the best crystals, (i) 
The evidence in the CMR context is overwhelming. Ex- 
periments and theory provide a picture where competing 
ferromagnetic (FM) and charge-ordered (CO) states form 
microscopic and/or mesoscopic coexisting clusters ||, ||. 
Exciting recent experiments M identified features re- 
ferred to as a "quantum critical point" (QCP) - de- 
fined as the drastic reduction of ordering temperatures 
near the zero temperature (T=0) transition between or- 
dered states - by modifying the A-site cation mean-radius 
(ta) by chemical substitution at fixed hole density (left 
inset of Fig. 1). The paramagnetic state in the QCP re- 
gion - where the Curie temperature Tc is the lowest ~ is 
crucial to understand CMR phenomenology, producing 
the largest CMR ratio ||i|, |, |l. (ii) In the HTS context, 
scanning tunneling microscopy (STM) studies of super- 
conducting (SC) Bi2212 revealed a complex surface with 
nm-size coexisting clusters Q . Underdoped cuprates also 
appear to be inhomogeneous pi. In addition, a "colos- 
sal" proximity effect (CPE) was reported on underdoped 
YBa2Cu306-(-x over large distances g]. 

In this paper, the competition between two ordered 
states in the presence of quenched disorder is investi- 
gated. These states are assumed sufficiently "different" 
that their low-T transition in the clean limit has first- 
order characteristics. The approach has similarities with 
the classical work of Imry and Ma [pi . From the general 
considerations, doped TMOs are here considered, with 
intrinsic disorder caused by chemical substitution. For 
Mn-oxides, a possible rationalization of the CMR effect 
is discussed, with predictions including a scale T* for 
cluster formation - the analog of the Griffiths tempera- 



ture [|lO[ but in the regime of competing orders. For un- 
derdoped Cu-oxides, a similar inhomogeneous picture is 
proposed. The calculations are mainly carried out using 
a two dimensional (2D) toy model of Ising spins, but simi- 
lar data in three dimensions (3D) and for the one-orbital 
manganite model have also been gathered. Then, our 
conclusions appear valid for a variety of models with com- 
peting orders. The actual Hamiltonian employed here, 
defined on a square/cubic lattice (spacing a=l) and with 
Ising variables, is H= Ji ^(ij) S^S^j + h E(im> S^SV 
+ J4X](in) S^iS^n, in a standard notation. Sites (ij) are 
at distance 1 (usual nearest-neighbors), (im) at distance 
■\/2, and (in) at distance -\/5. The three couplings are an- 
tiferromagnetic (AF). More than one coupling is needed 
to generate two competing T=0 states, and Ji and J2 
are the natural ones. However, the clean-limit first-order 
transition between those states was found to be more 
robust if a small J4~0.2Ji couphng is added. The result- 
ing competing states Oi and O2 are an AF state for low 
J2/J1, and a "coUinear" AF state with rows (or columns) 
of spins up and down for large J2/J1 [0. The main fea- 
tures of the toy model phase diagram are common to a 
variety of models with competing tendencies. 

The toy model phase diagram, without disorder, is 
shown in Fig. 1, and it has the expected shape: the or- 
dering temperatures decrease and meet at the clean-limit 
critical coupling J2c=0.7Ji, and the low-T transition was 
found to be clearly first-order. The most interesting re- 
sult in Fig. 1 is the influence of disorder on the clean-limit 
diagram. The first-order transitions become continuous 
with sufficiently large disorder, in agreement with pre- 
vious work P4 I. Critical temperatures far from J2C are 
not affected much by the disorder strengths considered. 
However, a drastic reduction is observed near J2c- In fact, 
the Monte Carlo (MC) results suggest that the obtained 
phase diagram is similar to the insets Fig. 1 for Mn- and 
Cu-oxides. With increasing disorder strength W, either a 
first-order line separating the competing phases still sur- 
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FIG. 1: Phase diagram of the 2D J1-J2-J4 toy model used to 
analyze the competition between ordered states. Results on 
32^ lattices (not shown) lead to a qualitatively similar phase 
diagram. Ji=l is the scale, and J4 is fixed to 0.2. Disor- 
der is incorporated such that J2(ij) at a link joining sites i 
and j is uniformly distributed between J2-W/2 and J2-I-W/2. 
Blue (red) curve corresponds to W=1.5 (0.75). Dashed black 
lines are the result without disorder W=0, and T* denotes 
the clean-limit transition. A Metropolis algorithm was used, 
on up to 512 lattices, calculating (i) the largest ordered clus- 
ter size [J. Hoshen et al, Phys. Rev. B14, 3438 (1976)] and 
(ii) the order parameters for AF and coUinear phases with a 
spin structure factor maximized at momenta (tt, n) and (vr, 0)- 
(0, 7r), respectively. Ten or more realizations of disorder were 
used, found to be sufficient for large systems. Insets are the 
phase diagrams of Mn-oxides in the FM-CO competition re- 
gion pj , and of the single-layer Cu-oxide in standard notation 
[J. B. Torrance et al, Phys. Rev. B40, 8872 (1989)]. Points 
1-5 are explained in Fig. 2. 



FIG. 2: AF order parameter vs T for the toy model at fixed 
W=1.5, using several lattice sizes with periodic boundary con- 
ditions, (a) corresponds to J2=0.65 and (b) to J2=0.69. Note 
in the latter the order parameter rapid suppression as the size 
grows, (c) Typical spin configurations representative of domi- 
nant 2D states. Shown are averages over 10 measurements, in 
about 100 MC sweeps to avoid correlations, after thermalizing 
with thousands of sweeps. Very similar results were obtained 
in 3D simulations. (1,2,3) are at J2=0.68, and T=2.00, 1.45 
(near the resistance peak, see Fig. 3), and 0.80, respectively 
(see Fig. 1). Green regions have coUinear order, while red and 
blue indicate Neel and "anti-Neel" order. The last two differ 
in the staggered order parameter sign, i.e. they intuitively 
a-re TiTi ••• a-nd itit .... The white does not have a dominant 
order after the MC sweeps considered here. Green/red/blue 
pale regions have weak order. (3', 4, 5) corresponds to T=0.8 
and J2=0.68, 0.70 and 0.73, respectively (see Fig. 1), and the 
Neel and anti-Neel states are here given the same color (red), 
while green remains coUinear. (d) Typical spin configuration 
at staggered field Hs=0.01, J2=0.68, W=1.5, and T=1.45. 



vives at J2C (red points), as in manganites, or a disordered 
region of finite J2 width opens at T=0 (blue points), as 
in single-layer cuprates. Note that the ordering tem- 
peratures exactly meet at T=0 for only one fine-tuned 
W. However, the overall shape of the phase diagram re- 
sembles QCP behavior in a robust range of W. For this 
reason, our results are qualitatively described as induc- 
ing "quantum-critical-point like" behavior, not a rigorous 
expression but hopefully descriptive enough to be useful. 

Sufficiently strong quenched disorder will smooth first- 
order transitions |1^. At J2C, this should occur at in- 
finitesimal W in 2D M, while a finite W is needed in 
3D. Then, at J2C and with finite temperature, a param- 
agnetic state must be generated with growing disorder. 
Note also that our toy model is classical, but QCP-like 
behavior is nevertheless generated p3| . 

Figure 1 is the result of a systematic computational 
effort. As example, in Fig. 2a-b, the AF order parame- 
ter vs T is shown for J2 values outside and inside the 
coupling range where a T=0 disordered regime is ob- 
tained. For J2=0.69 note the order-parameter cancel- 



lation with increasing size (this coupling is not critical at 
Tt^O). Representative spin configurations are in Fig. 2c. 
Keeping J2 constant and varying T, three regimes are 
observed: (1) A high-T regime, where the system is dis- 
ordered after MC time averaging. (2) An intermediate 
region Toi<T<T* with preformed clusters, but with un- 
correlated order parameters giving a globally paramag- 
netic state, similar to the Griffiths phase. (3) A low-T 
regime where the clusters from (2) grow in size, although 
the disorder is uncorrelated from link to link, and perco- 
late upon cooling. Note that clusters with different signs 
for the order parameter are separated by thin regions of 
the competing phase, providing a possible mechanism for 
stabilizing domain walls. Considering now fixed low-T 
but changing J2, configurations (3', 4, 5) are obtained. In 
this case, just the two ordered phases are in competition 
(no white regions), and the transition between phases 
appears percolative-like. 

The main features of the results Figs. 1-2 - shape, 
clustered structure, the QCP-like behavior - are believed 
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FIG. 3: (a) Proposed state for Mn-oxides in the CMR regime, 
(b) Resistance of the toy model after the equivalence to man- 
ganite states is used (see text), using a 256^ lattice, W=1.5, 
at the couplings indicated. The calculation is carried out by 
transforming a spin configuration into a resistor network, with 
nodes centered at plaquettes (# nodes = 1/4 # sites) and 
resistors between them. The values of the conductances of 
these resistors were established using BB=RR=1.0, BR=0, 
WW=0.3, GG=BG=RG=WG=0, WB=WR=0.5, where B, 
R, G, and W stand for blue, red, green and white regions (see 
Fig. 2c), AA'= a means that the resistor between plaquettes 
in the A (=B,R,G,W) and A' states has value a (arbitrary 
units), and AA'=A'A. Other values for WW and WB lead to 
similar results, and BB defines the scale. Note that the con- 
ductivity should be spin dependent, and a BR link (when an 
electron moves from a spin-up to a spin-down region) has zero 
conductance. The algorithm used to obtain the total conduc- 
tance is exact |l4]. (c) Resistance (arbitrary units) vs T, at 
external fields Hs indicated, using a 256^ lattice, J2=0.68, and 
W=1.5. (d) Tq results for the one-orbital manganite model 
using 8^ and IG'^ clusters, density x=0.5, infinite Hund cou- 
pling, and hopping t=l. With disorder, the Jaf couplings are 
randomly distributed between Jaf-W and Jaf+W (W indi- 
cated). In practice, Tq was defined when the spin correlations 
at distance \/2 dropped below 40% of the maximum value. 



to be qualitatively general. In fact, simulations of one- 
orbital models (below) and other models studied in this 
effort give a similar phase diagram. Of course, the anal- 
ogy should not be taken too far, e.g. critical exponents 
may not be universal since the Ising model underlying 
symmetries are quite different from those of realistic sys- 
tems. However, it is worth investigating the consequences 
of the general phase diagram Fig. 1 for materials where 
two states strongly compete, such as in Mn-oxides. In 
this context, if a simulation of a realistic model with FM 
and AF phases on a huge lattice were possible, FM and 
AF clusters analogous of Fig. 2 would be found. Then, 
a reasonable way to bypass that (currently impossible) 
computational effort is to simply translate Fig. 2 into 
manganite language. This is a speculation, but hopefully 
the essence of the problem is preserved by the procedure. 
The proposed translation links order Oi with ferromag- 



netism, with order parameters pointing in different di- 
rections for different clusters, while O2 corresponds to 
charge-ordering. Translating Fig. 2c(2) into Mn-oxide 
language leads schematically to Fig. 3a, our proposed 
CMR state. The preformed FM clusters have uncorre- 
lated moment orientations, and zero global magnetiza- 
tion. Note also that the "depth" of the QCP-like feature 
is not universal, it depends on the disorder strength. 

To test the relevance of Fig. 3a to CMR manganites, 
a resistor network calculation was set up. Translating to 
Mn-oxide language, as explained before, the MC gener- 
ated configurations were mapped into a resistance grid 
(see caption of Fig. 3b). For up (down) spins, the blue 
regions of Fig. 2c - analog of positive magnetization FM 
clusters - have high (low) conductivity, the red regions 
have low (high) conductivity, and the green regions are 
insulating. The Kirchoff equations were solved exactly, 
leading to the results Fig. 3b. In agreement with in- 
tuition, the non-percolated state Fig. 2c(2) has a large 
resistance for both spins up and down, while the perco- 
lated low-T or disordered high-T states have far better 
conductance. Note that the resistance peak intensity in- 
creases as the ordering temperature is reduced varying 
J2, analog of (r^i), closer to the QCP-like regime. 

The rotation of the large moments of the preformed 
FM clusters (Fig. 3a) may occur with small magnetic 
fields. These effects are mimicked in the toy model us- 
ing a staggered external field Hg which favors Oi clusters 
with order parameter Ms>0 (blue, Fig. 2c) to the detri- 
ment of Oi clusters with -Mg (red, Fig. 2c), or O2 regions. 
Figure 2d confirms the rapid generation of positive Oi or- 
der in the region Tqi <T<T* with tiny fields O.OlJi. This 
severely affects transport (Fig. 3c), i.e. a modest field 
transforms the intermediate T cluster state into a fairly 
uniform state with robust conductance. Results Fig. 3c - 
the main results of this paper - are similar to those found 
in Mn-oxides, and a huge MR ratio [R(0)-R(IIs)]/R(IIs) 
of-4xl03 % was obtained at Hs^lQ-^ |i|. 

In addition, there are already experimental indications 
in Mn-oxides for the existence of a temperature scale T* 
for uncorrelated cluster formation [|l6| , which should be 
ubiquitous in low-bandwidth manganites ||l3l ■ 

The phase diagram Fig. 1 is representative of more re- 
alistic models. Fig. 3d contains MC results for the one- 
orbital model p8[ , in which the coupling Jaf between 
localized spins is varied to induce a metal-insulator tran- 
sition 0. Without disorder, the T'-^^O transition is known 
to be first-order between FM and AF states, the latter 
with rows or columns of spins up and down f^. The 
"characteristic" ordering temperatures TJ. , at which spin 
correlations become robust upon cooling, are shown vs. 
Jaf- Note the similarity with Fig. 1. 

The results in Figs. 1-2 can also be adapted to 
cuprates. The high-Tc phase diagram (inset Fig. 1) 
shows a suppression of AF and SC order in a region usu- 
ally labeled "spin-glass" , whose origin is unclear. Consid- 
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FIG. 4: (a) Conjectured HTS phase diagram. Black lines 
should be the actual phase boundaries without disorder. The 
shaded region is conjectured to have metallic (SC) and insu- 
lating (AF) coexisting regions in the real materials, (b) Re- 
sistivity p„h ys T, from a random-resistor network calculation 
as in Ref.pd], where details can be found. A 50x50 cluster 
was used, with pab for insulating (optimal doping) fraction 
p=1.0 (0.0) taken from LSCO x=0.04 (0.15) data [H. Takagi 
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ering these diagrams together with the CPE results [pj, 
it is conjectured that the very underdoped cuprate state 
may not be homogeneous but results from a SC vs. doped 
AF competition after disorder is considered. Inhomo- 
geneities (clusters) should be present even within ordered 
phases [Fig. 2c(3)]. Stripe states are candidates for the 
doped AF state ^. 

The proposed clean-limit phase diagram is in Fig. 4a, 
with a vertical first-order transition line, as cuprates have 
upon electron doping |20(| , and heavy fermions with vary- 
ing pressure. The shaded region may contain a mixture 
of stripe-like and preformed SC islands |2^. Due to the 
general character of the discussion of Figs. 1-2, colossal 
effects should he ubiquitous when ordered phases compete, 
and they could appear in cuprates as well. CPE [gj could 
be a manifestation, with preformed SC clusters percolat- 
ing under the influence of nearby SC materials. To fur- 
ther check this hypothesis. Fig. 4d contains results of a 
phenomenological random-resistor calculation of resistiv- 
ity vs T p2[ in rough agreement with experiments p3|. 

Summarizing, the results presented in this paper |24| 
suggest that "colossal" effects in TMO's could originate 
in intrinsic inhomogeneities. These large effects may be 
more general than previously anticipated. In our stud- 
ies, the analog of the classical Griffiths regime - usually 
associated with weak effects - is here much more robust, 
strongly susceptible to external fields. 
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